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ABSTRACT 

We develop the Callan- Witten picture for monopole catalysed skyrmion 
decay in order to analyse the corresponding cosmic string scenario. We discover 
that cosmic strings (both ordinary and superconducting) can catalyse proton 
decay, but that this catalysis only occurs on the scale of the core of the string. 
In order to do this we have to develop a vortex model for the superconducting 
string. We also give an argument for the difference in the enhancement factors 
for monopoles and strings. 
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Several years ago Callan [l] and Rubakov [2] showed how it was possible for a grand 
unified monopole to catalyse proton decay. They also showed that the catalysis would occur 
with an enhanced cross-section: the inverse square of the proton rather than the grand 
unified mass scale. That a grand unified topological defect can catalyse baryon decay 
is not surprising. In such objects, the symmetry of the grand unified group (e.g. SU(5) 
or SO(IO)) is unbroken in the core of the defect, thus allowing baryon non- conserving 
reactions to occur. What is somewhat surprising is the enhancement of the cross-section. 
Naively we might expect that only a ‘direct hit’ of the proton on the monopole core would 
allow such a catalysis to proceed, however, due to the magnetic moment of the fermion 
there is an attractive force binding it to the monopole - thus leading to the enhancement 
factor. 

The cosmological implications of this result are immediate. If topological defects exist, 
then they can catalyse baryon decay. If they catalyse baryon decay at a sufficient rate then 
the effects should be observable. The lack of observational data to support proton decay 
at such rates then places a constraint on the number density of such defects. Thus the 
enhancement factor in the case of monopoles actually places more stringent bounds on the 
monopole flux [3,4] than conventional methods [5]. 

Cosmic strings are also topological defects of grand unified theories. Thus, like 
monopoles, they can catalyse baryon decay. However, it has recently been shown that 
for both cosmic [6] and superconducting cosmic strings [7], there is no enhancement of the 
cross-section; the catalysis cross-section is just given by the grand unified scale. Since the 
density of strings is very high at early times in the evolution of the universe, it is still 
possible that string catalysis could have cosmological implications despite the smallness of 
the cross-section. For example, string catalysis could place constraints on baryogenesis [8]. 

It is therefore important that monopole and string catalysis are fully understood. The 
quark-monopole(string) scattering is fairly well understood, but it describes a high energy 
process. In low energy processes (such as decay) we should consider a more appropriate 
hadron model. In the cosmological context, where we expect fairly low velocities, it is 
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important that a more appropriate model is developed. One such model is the Skyrme 
model [9], which has the advantage over hadronic models in that it is easy to couple to 
strings or monopoles. A more realistic low-energy hadron model, such as the skyrmion, 
might yield some additional insight into the physical decay process. Callan and Witten [10] 
were the first to examine such a setup. They considered a skyrmion model for the hadron, 
a Dirac monopole, and ex ami ned the processes occurring during the skyrmion-monopole 
interaction. 

We use this as a starting point for our discussion. In section 2 we give a brief resume 
of the Skyrme model. We review and develop the Callan- Witten [10] argument in section 3 
using the Wu-Yang [11] picture of a monopole. Such a picture removes the problem of the 
distinction between the physical singularity of the electromagnetic fields at the monopole 
core and the Dirac string singularity, which is a gauge artefact. Our argument also clarifies 
aspects of the topological unwinding of the skyrmion on the monopole. In the fourth section 
we examine the scattering of a skyrmion off a cosmic string. We first use the wire model 
for the string in order to mimic the Dirac model for the monopole, however such a picture 
does not permit baryon decay. We are therefore forced to consider a vortex model for the 
string in order to obtain catalysis in the string core. In section 5 we consider the analogous 
process for a superconducting string [13]. First we use the wire model, but despite there 
being a long range force in this case, we again show that such a picture does not result in 
baryon decay. We then develop a vortex model for the superconducting string, by solving 
the string equations of motion near the core. In the vortex model we obtain catalysis in 
the string core. Our analysis gives a heuristic explanation of the enhancement factor with 
monopoles, which we explain in the concluding section. 
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2. Resume of the Skyrme model. 

In this section we review the Skyrme model of the nucleon [9] . This is a sigma model 
with stable soliton solutions otherwise known as skyrmions. (For a good review see e.g. 
Balachandran [14].) The lagrangian is 

Csk = i/5 Trd„Vd»V-' + jijTtJ V-'d,.V,V-'drU]' , (2.1) 

where a is a constant of order one. Here U in an SU(2) matrix related to the pion fields 
via 

U = exp{^- r.7r} , (2.2) 

where r = (ti,T 2 ,T 3 ) are the three generators of SU(2). An alternative representation of 
U, which will be useful in what follows, is one which makes explicit reference to the S 3 
topology of 517(2): 

U = cos + i sin ^—^r.7r . . (2.3) 

fir fir 

U is now represented as a unit four- vector in terms of the basis {/, ,r}, where J s is the 
identity matrix. 

The lagrangian (2.1) admits soliton solutions, i.e. localised stable finite energy field 
configurations. A topological picture of this can be drawn in the following way. Finiteness 
of the energy requires that U(x) — f const, as |z| -* oo. We can therefore think of a 
soliton field configuration as a map from compactified three-space (IR 3 U{oo} = S 3 ) to the 
three-sphere of SU(2), 

Usol (*) : Slurs ^«/(>) - ( 2 ' 4 ) 

Such maps may be classified according to the homotopy equivalence class to which they be- 
long. Members of the same equivalence class are related to each other by a 
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continuous deformation and represent translated or excited states of the same soliton. 
Since Il 3 (S 3 ) = 7L, we may conclude that soliton field configurations are labelled uniquely 
by an integer value, Nb (the baryon number), which is the degree of the map (2.4). In a 
dynamical theory, the continuity of the fields implies that Nb is a continuous function of 
time and hence constant. 

The baryon number may be alternately represented as the charge associated with the 
conserved current 


B* = -^e^Titu-'dvUU-'dMU-'dvV) 

5 247 r 2 \ / 

N b = J B°d 3 x . 


(2.5a) 

(2.56) 


The standard nucleon field configuration ( Nb = 1) is usually represented by the 
Hedgehog configuration: 

tfjv(*) = exp [tF(r)i .t] , (2.6) 

where -F(O) = w and F(oo) = 0. In practise F(r) differs significantly from zero only within 
the core of the skyrmion, i.e. a distance ~ (a/*) -1 . 

In the presence of electromagnetism, the preceding discussion must be modified to 
allow for the nucleon charge and magnetic moment interaction. We must generalise Csk 
to be invariant under 


U-*U+ iea[Q , U] 


(2.7a) 

(2.76) 


where Q is the quark charge matrix: 


-(! A) 


( 2 . 8 ) 


(For simplicity we are considering the case of two quark flavours only.) 


4 



ORIGINAL PASS VS 

OF FOOR QUALITY 

Taking into account QCD anomalies, the lagrangian and baryon current become, [15], 
C = ±f} TiV^VW-' + 

+ -^—e^ApTrQldvUU-'dpUU-'drUU- 1 

l07T 2 

+ u- 1 d u uu- 1 d p uu~ 1 d IT u) 

+ ^-rt^^A^ApTt^daUU- 1 + Q 2 U~ 1 d a U 
87 r 2 L 

+ iQd'UQU- 1 - \QUQd v U~^ (2.9o) 

B p = B p + -^c* ,/p,r d v [A p TiQ{U- 1 d l ,U + drUV- 1 )} , (2.96) 

o7T^ 

where V^U = d^XJ — ieA^lQ, U] is the covariant derivative of U. 

Note that (2.7b) is infinitesimal in form. We will need the non-infinitesimal version: 

U e iea(x)Q Ue -iea{z)Q 


Noting that Q = |J, + jr 3 , we see that this reduces to 


( 2 . 10 ) 


Note that the generalisation of the baryon current (2.9b) now contains a term depen- 
dent on Ap. This appears via a term which is a divergence, thus provided there are no 
singularities in A and that surface terms vanish, the baryon number is still integral. In 
terms of the topological picture presented previously, provided there are no singularities, 
U SOL (x) is still a map from S 3 — » 5 3 and thus the classification of maps into equivalence 
classes labelled by baryon number still holds. 
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3. Monopole catalysed skyrmion decay. 


In this section we review the argument of Callan and Witten [10] for skyrmion decay. 
Our presentation is from a different viewpoint in order to facilitate the transition to the 
string picture. For simplicity let us now consider the Dirac monopole. In its original 
formulation this has a gauge potential given by* 

A<p = 5(1 - cos 6) . (3.1) 

In the ensuing discussion, we will assume that eg = 1/2, in accordance with the Dirac 
quantisation condition. This is singular on the line 9 = 1 r, however the electromagnetic 
flux, 

Fg <t> = g sin 0, (3.2) 

is finite everywhere except at r = 0. Thus the singularity of A^ on 9 = ir is a gauge artefact, 
the Dirac string. It arises because we are trying to express the electromagnetic field tensor 
as the exact differential of a covector gauge field on IR 3 — {0}. Since the gauge field now 
does contain singularities, the interpretation of the divergence in the second term of (2.9b) 
is unclear - is the baryon number still integral? The mathematical/topological picture 
also becomes clouded. If one removes the semi-infinite singular line { 0 = 7r} from IR 3 , the 
soliton field configurations are now maps from IR 3 — = 7r} (which is contractible) into 5 3 . 
These are all topologically trivial. Clearly this is too naive, for the field U must satisy some 
continuity property near 9 = ir in order for it to contain no singularities. However, since 
U is coupled to the gauge field A ft, we must be careful in specifying boundary conditions 
along 0 = ir. 

In order to make these intricacies more transparent, we will take an approach to the 
Dirac monopole which avoids Dirac strings — that due to Wu and Yang [11]. Briefly, the 
singularity in A M is exactly analogous to the coordinate singularity we obtain if we try 
to cover S 2 with only one coordinate patch. Just as this singularity can be removed by 

t Note that here, and throughout, we use a coordinate basis for our components and 


not cm orthonormal one. 
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choosing two patches, the singularity in A ^ can be removed if one chooses instead two 
coordinate patches for IR 3 — { 0 }, each with an associated A^, relating the two different 
‘branches’ of A^ by a gauge tranformation on the overlap. 

Two convenient patches are 


(0 

0 < 0 < 7 r — 6 ; r >0 

(3.3a) 

(*) 

6 <0 <ir ; r> 0 , 

(3.36) 


with 


= $(1 - cos 0 ) 

(3.4a) 

Arf = -< 7(1 + cos 0 ). 

(3.46) 


These are related by the non-trivial gauge transformation 

A t n = A l f i — 2gd tl <f> (3.5) 

on the overlap. The Dirac string picture is the limit as 6 — » 0 of one coordinate patch. 
We see that this picture now has no coordinate singularities. The gauge field on each 
coordinate patch is perfectly regular. The electromagnetic flux is g sin 0 on each patch and 
is independent of the branch of A ^ chosen on the overlap. This is what we would expect 
for a physical quantity. 

We now want to include the SU( 2 ) field, U, in this picture. The nucleon field is 
non-singular, therefore it must be well defined everywhere. Since the U -field is coupled to 
the gauge field, the presence of the two branches of Af, indicates that we must define a 
separate field configuration on each chart. These will then be related in the overlap by a 
non-trivial transformation induced by the gauge transformation (3.5) on A M . From (2.10) 
we conclude that this is 

V , = e ”*^ T * / 2 17 x e*^ T * I 2 . (3.6) 

We now have a perfectly consistent, singularity free picture of the nucleon on the back- 
ground field of the monopole. 

Having removed the singularity problem, we see that once again the SU( 2 ) field config- 
uration is a map from compactified physical space into the SU( 2 ) three-sphere. However, 
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whereas before we could regard IR 3 as the union of the two coordinate patches (l) and 
(t ) with the trivial gauge transformation on U, here we have a non-trivia! transformation 
for U in the overlap. Thus although we can classify the field configurations in each case 
according to homotopy equivalence, there is no reason to assume that in each case these 
classes will be the same. Indeed, if we write U in the form (2.2), the effect of the gauge 
transformation (2.10) is to rotate the vector w by an angle <f> around the 3-axis. This will 
have a twisting effect on the i ?\ , w 2 components. Thus the presence of the monopole gauge 
field shuffles the members of the equivalence classes of (2.4). This ‘shuffling* is crucial to 
the physical description which follows. 

For convenience and comparison with the cosmic string case we will subdivide the 
argument into three steps: firstly we show that there is a suppression of the wave function 
of a charged particle near the monopole core; secondly, that due to the shuffling of the 
baryon equivalence classes it is possible for a pure 7r° field configuration to carry baryon 
number; and thirdly, that under the presence of suitable boundary conditions there can be 
a non-zero radial baryon flux — the skyrmion unwinds. 

i) Suppression of charged particles near core. 

Solving the Klein- Gordon equation, 

(V„V" - e'ADv = 0, (3.7) 

for the charged pions in the presence of a magnetic monopole shows that charged pions 
can no longer have zero angular momentum. Instead of the usual spherical harmonics, we 
must now use generalised monopole harmonics [12] for the wave function in the monopole 
background. These have a dependence on the charge of the particle. The lowest harmonic 
for the patch (1) is <p oc y/l + cos 6. Substituting this into (3.7) gives for the radial part 
of the wave equation: 

1 _ 2o 1 

_ 5rr drip = —y, . 

This implies an asymptotic behaviour of ip cx r » as r — » 0. Thus the wave function of 
the charged pions is suppressed near the core. Note however that for uncharged particles 
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no such suppression occurs. Thus in order for the nucleon to approach the monopole core, 
it must be able to deform into a pure 7r° field configuration. In order for this process to 
be possible, the 7T° field configuration must be able to carry baryon number. 

ii) Baryon number of a 7r° configuration. 

For simplicity we consider a purely radial 7r° field configuration 

Uk = exp |*/t 3 J 

= cos/ + tsin/r 3 , (3.8) 

where /( 0) = 0 and /( oo) = 27 r. It is then straightforward to calculate the baryon number 
from (2.9b): 

Btt = ^rHv-'WV-'Wu-'arv) 

+ + drUU- 1 )}, ( 3 . 9 ) 

and the expression for the derivative of Ujc : 

dnU K = / )M (-sin/-fzcos/r 3 ) • 

= if^T 3 U K . (3.10) 

Clearly the first term vanishes, since / depends only on r and possibly t. The 

second term gives 

= + vv-') 

= "j 

r 

87r 2 r 2 

=*■ ■ Wb= ^M.” = 1 ’ < 3 - u ) 

thus demonstrating that it is possible for a pure 7r° configuration to carry baryon number 
in the presence of the monopole [10], The configuration (3.8) is called the radial kink [10]. 
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iii ) Unwinding of the radial kink. 

Now let us consider the radial baryon flux of the kink. B r vanishes, since / depends 
only on r and t. Therefore 


B r = ^ t TV(, °d v [A p TiQ(U^ d„U K + d a U K U^)) 

= £?h‘ > ‘ A * TTe>(V « i ' K + 

m f_ 

87r 2 r 2 


(3.12) 


Thus we see that the radial flux of baryon number into the monopole core is 2*^' Whether 
or not /(0,t) can be non-zero depends on the boundary conditions at the monopole core. 
In the case of a grand unified monopole formed during an SU(5) or SO(IO) phase transition 
for example, it is possible for baryon non- conserving boundary conditions to be placed, 
and hence for /(0,<) ^ 0 [10]. 

Thus monopoles can catalyse skyrmion decay. We will now give a brief summary of 
the essential dynamical process. 

Let us suppose that a standard nucleon field configuration is approaching the monopole 
and that the monopole will encounter some region in which the field configuration is non- 
trivial. Because of the suppression of charged particle wave functions near the monopole 
core, the nucleon field configuration will be forced into a radial kink. After passing the 
monopole the nucleon will regain a more conventional field configuration profile. 

In passing through the core, the skyrmion must make the transition from one coordi- 
nate patch to another. Thus in the absence of baryon non- conserving boundary conditions, 
the field configuration must ‘twist’ as it passes over from one coordinate patch to another. 
This is in order to preserve baryon number. If, on the other hand, the boundary conditions 
do admit baryon decay, the skyrmion need not twist, but merely keeps its original profile, 
which in the new coordinate patch has baryon number zero. 

This is the picture for monopole catalysis of skyrmion decay. 
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4. String catalysis of skyrmion decay. 


We now turn to the case of a skyrmion interacting with a cosmic string. At first 
sight, we might expect some similarities with the monopole case, since the monopole has a 
semi-infinite Dirac string singularity, and we have an infinite string. However, this would 
be misleading; the Dirac string is a gauge singularity and can easily be removed by a more 
suitable description in terms of coordinate patches. 

This is the crucial difference from the monopole setup. For a non-singular approach 
to the gauge field outside a monopole we needed to define two branches of the gauge field 
on two different coordinate patches, related by a non-trivial gauge transformation on the 
overlap. The cosmic string however, has a perfectly well defined gauge field without invok- 
ing coordinate patches. Thus the gauge field for a cosmic string exhibits no singularities, 
the additional term in (2.9b) is once more a total divergence, and baryon number is un- 
changed. Alternatively, if there are no gauge singularities, the equivalence classes of the 
soliton maps (2.4) are unchanged. 

We will consider the equivalent of steps i) to Hi) for the case of cosmic strings firstly 
by using a wire model for the string, and then by a more realistic vortex model. 

I. WIRE MODEL 

In grand unified models the string width is of the order of M ~* , where M is the grand 
unified mass. Thus, as a first approximation we take the string as a wire singularity [16] 
on the symmetry axis. In this case the string energy-momentum tensor is represented by 
the distributional form 



Away from this singularity the gauge field is given by 


Afj. — 



in cylindrical polar coordinates {p,i?,z}. 


(4.1) 
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In cylindrical polar coordinates, the static Klein-Gordon equation (3.7) reduces to 

(V„V“ - t>Al)v = - + a l + ^( a l - «’•*})]*> 

= 0. (4.2) 

Here, rather like the monopole case, <p picks up extra “angular momentum” around the 
2 -axis due to the presence of a non-zero A#. For the wire model, (4.2) gives 

f - a pP a, + a] + i(«J - i)l v = o , ( 4 . 3 ) 

Ip p* J 

which implies that the radial part of the wave equation for the lowest angular momentum 
eigenstate is 

pd p pd p ip(p) = <p{p). (4.4) 

From this we deduce that <p must tend to zero as least as quickly as p near p = 0. Therefore, 
as in the monopole case, the wave functions of charged particles are suppressed near the 
core of the string, but those of uncharged particles need not be. 

Unfortunately, it is now impossible for a radial kink to carry baryon number as a quick 
glance at (3.9,10) shows. B £ is zero as before, hence 

B° = ■^z€ 0vpr d v [ApTTQ(U~ 1 U' + U'U- 1 )). (4.5) 

OTV 1 

Translating (4.1) into spherical polars, we obtain 

Ap = > ( 4 * 6 ) 

t 

which is a constant. We can now see immediately that B° = 0. 

We have run into a problem here. Taking the wire approximation for a cosmic string 
leads to a suppression of the charged pion fields near the string. However, since a radial 
kink cannot carry baryon number on this case, we cannot have a deformation of the nucleon 
fields that would allow the skyrmion to approach the string core. Hence in the wire model 
of cosmic strings we do not get catalysis. Perhaps this problem is a result of approximating 
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the string core by a line. In order to be more physically realistic, we will consider a vortex 
model for the string — the Nielsen- Olesen vortex. 

II. THE NIELSEN-OLESEN VORTEX 

To illustrate the salient features of skyrmion catalysis by cosmic strings it is only 
necessary to consider an abelian theory. Thus we consider the Nielsen- Olesen vortex [17]. 
We discuss this in detail to facilitate the extension to the superconducting case. The 
Nielsen- Olesen string is a vortex solution to the lagrangian 

£[*, A„] = - 1 F ll ,F‘“‘ - jWV - (4.7) 

where + ieA^ is the usual gauge covariant derivative, and F^ v the field strength 

associated with A M . However, we shall choose to express the field content in a slightly 
different manner, and one in which the physical degrees of freedom are made more manifest. 


We define the (real) fields X, x and by 

<£(*“) = rjX(x a )e ix ^ a ^ (4.8a) 

M*’) = ; fr(* a ) - V„x(**)] . . (4.81) 

c 

In terms of these new variables, the lagrangian becomes 

L = rfV^XV^X + t fX'P^P* - 1~F^F* V - X 2 - l) 2 (4.9) 

4e 2 4 

where is the field strength associated with P ^ . The equations of motion are 

SS 

— = 2r) 2 V^X - 2r l 2 P li P»X + \rfX (X 2 - 1) = 0 (4.10a) 

- 2X 2 i) 2 P v = 0 . (4.106) 


We can see that a vacuum state is characterised by X — 1. However, this is not the 
only stable ground state which solves the equations of motion. Nielsen and Olesen showed 
that there exists a non-trivial stable ground state solution to the above equations of motion 
which has a vortex-like structure. 
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The Nielsen and Olesen vortex solution corresponds to an infinite, straight static string 
aligned with the z-axis. In this case, we can choose a gauge in which 

<f> = nX{p)e ie ; P" = P(p)V"0 (4.11) 

in cylindrical polar coordinates. This string has winding number one. The equations for 
X and Pft greatly simplify to 

-X" - y + ^ + —-Xtx 2 - 1) = 0 (4.12a) 

-P" + ^+ 2 e 2 rj 2 X 2 P = 0 . (4.126) 

These equations do not have any known analytic solutions, but asymptotic forms may 
be derived. It is not difficult to see that these are: 

X oc p ; P = l + 0(p 2 ) as p -> 0, (4.13a) 

X a 1 - p -i/ 2 c -^p . p a ^fpe-^P as p — » oo. (4.136) 

These will be sufficient for our purposes. 

We now wish to use this Nielsen- Olesen vortex as a model for the cosmic string. Thus 
instead of using the wire form for the gauge field, (4.1), we use the expression for the 
Nielsen- Olesen gauge field, which is given by (4.8) and (4.11): 

K = \(P(p) - . (4.14) 

This modifies the Klein-Gordon wave equation. The radial equation for the lowest angular 
momentum eigenstate, corresponding to (4.4), becomes 

pd p pd p <p(p) = (P(p) - l)V(p) . (4.15) 

From (4.13a), P(p) = 1 + 0(p 2 ) near p = 0, hence 

pd p pd p <p(p) = 0(p 4 )yj(p), 
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allowing y? ~ const, as p — * 0. Thus, on the scale of the core of the string, we need 
not have total suppression of charged particle wave functions. 

Writing the analogue of (4.6) for the vortex A ^ in spherical polar coordinates gives 

A m = ~(P(r sin 6) - 1)V ^ . (4.16) 

e 

As before, although slightly less trivially, substituting this into (4.5) shows the radial kink 
cannot carry baryon number. To see this, note that 

B° = ’' pr d u [A p TrQ(U- 1 U ' + U'U- 1 )) 

o7T^ 

= ***™x*r*) 

__ f dgP(r sin 6) 

47 r 2 r 2 sin 8 

=» N b = -^J o <fr/'[p(rsinfl)]* = 0. (4.17) 

However, this is no longer critical for we can have all three pion fields approaching the 
core. Once the skyrmion is in contact with the core of the string, where the grand unified 
symmetry is essentially restored, the possibility of decay arises. 

We will consider an unwinding process involving all three pion fields by making the 
simple Ansatz that the nucleon field configuration now depends on time: 


UN{z,t) = exp[iF(r,<)i.T]. ( 4 - 18 ) 

The calculation of the baryon current for this field configuration is somewhat involved, 
and we relegate the details to an appendix. 

The main result we need is the radial baryon current of the field configuration (4.18) 
which, from the appendix ((A15)), is given by 


B r 


F 

47 r 2 r 2 


[P(cos 2F — 1) + rP' —r—jr] . 

sin? J 


(4.19) 
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Integrating this over a sphere of radius r gives 

f B r d 2 x = —^~ [ dd sin 6 [(cos 2F -l)P(r sin 8) + r^r~lrP'(T sin 0)] 
Js * 27T J 0 sin a 

= / d0 [(P cos 0),$ + P cos 2P sin 0] 

2?r y 0 

J dd sin flP(r sin0). 


P Pcos 2P 

7T 271 


(4.20) 


For small r, P(r sin 0) = 1 + 0(r 2 ) implies 


J. 


B r d 2 x ~ — (1 — cos2P) 

5» 7T 


(4.21) 


and hence that the flux of baryon number into the string core is — F(1 — cos2P)/tt. 

Thus in the presence of baryon non-conserving boundary conditions, such as we would 
expect in the string core where the grand unified symmetry is unbroken, the skyrmion can 
unwind. Since F(0) = ir and F(oo) = 0 for the standard nucleon field configuration (2.6), 
we expect that for an unwinding process F changes from it to 0 at the core of the skyrmion. 
And indeed 


A Nb = J Bffdt 

~ J dtF(l — cos2P)/7r 
= — A [F — £ sin2F] 


7T 
= -1 


The residual field configuration is a topologically trivial excitation of the pion fields, and 
can therefore dissipate. 

Thus strings can catalyse skyrmion decay. The picture however relies fundamentally 
on taking a vortex model for the string, i.e. one in which the string has a finite thickness. 
A model of the string with infinitesimal thickness (the wire model) gives no catalysis. 
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5. Skyrmion catalysis by superconducting cosmic strings. 

In this section we consider the case of the skyrmion interacting with the superconduct- 
ing cosmic string. Unlike its Nielsen- Olesen cousin, this has a long-range electromagnetic 
gauge field, hence we might expect some differences with the previous analysis. After all, 
one of the differences between the monopole and the Nielsen- Olesen vortex was the absence 
of long range interactions in the latter setup. However this is not the case as we will now 
show. First we consider the wire model for the superconducting string, then we derive a 
more realistic vortex model by solving the equations of motion near the core. 

I. THE WIRE MODEL 

Similar to the cosmic string case discussed previously, we can try taking the super- 
conducting string to be a wire singularity on the symmetry axis. The long range electro- 
magnetic gauge field is 

A *(p) = 2“log(/>/p 0 ), O 5 - 1 ) 

where po is the radius of the string, and I is the current in the string. Imposing (5.1) for 
p > 0 gives a wire model for the superconducting string. 

Since we now have a long range electromagnetic field, we might expect some modifi- 
cations of the previous analysis. Consider first the Klein-Gordon equation. In cylindrical 
polar coordinates, the Klein-Gordon equation (3.7) reduces to 

(V,V<* - e ! 4V = - [- p S f pd„ + (82 - e*Al) + i«S] V 

= 0. (5.2) 

Thus, similar to the monopole and cosmic string cases, <p picks up extra “angular momen- 
tum” due to the presence of a non-zero A x . When we insert the form for A z from (5.1) 
into (5.2) there is no analytic solution for <p. However, it is possible to show that charged 
particle wave functions are suppressed near the wire, but those of uncharged particles need 
not be. 

In order to see if the radial kink can carry baryon number we need to consider the 
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corresponding form for A ^ in spherical polar s, 

A r = cos 6A x {t sin 6) 

Ae = — r sin 6 A z (r sin 8) (5-3) 

A# — 0 . 

As in the previous discussions, B £ is zero for the radial kink. From (3.9) and (3.10) we 
obtain 

= ~J~2 9 " A - Tt W k U 'x + V 'k u k )• < 5 ' 4 ) 

and we see once again that the baryon number of the radial kink (3.8) must be zero. This 
is because there is no <£- component or ^-dependence in A^. Thus the previous discussion 
given for the ordinary cosmic string also applies to the case of superconducting cosmic 
strings: since the charged fields cannot approach the string core, and since a radial kink 
cannot carry baryon number, the nucleon cannot approach the core and unwind. In order 
to get catalysis, we have to look at a thick vortex model. 

Incidentally, if we take a wire model of the U(l)xU(l) bosonic superconducting string, 
by setting the non-electromagnetic gauge field to the Nielsen- Olesen form (4.1), the be- 
haviour of <f> close to the wire is dominated by this field. The radial suppression of <f> 
becomes exactly the same as for non-superconducting strings. 

II. VORTEX MODEL FOR THE SUPERCONDUCTING STRING. 

In order to obtain catalysis it seems necessary to consider a vortex model for the 
superconducting string. To obtain such a model, we consider the U(l)xU(l)' model of 
Witten [13]. Previously, only numerical solutions have been found [18], so we give details 
of our solution. Our discussion follows closely that for the Nielsen- Olesen string discussed 
in the previous section. 

The lagrangian in this case is 


C = + D r <r'D»a - 


- 1 1 ) 2 + (M* - 


(5.5) 
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where <f> and <r are complex scalar fields; and / are coupling constants; 

Vn<f> = V n<i> + igC h 4> 

V^cr = Vft<r + ieA^c, 


Cft and An being abelian gauge fields carrying charges of g and e respectively, with 
and Fp V being the corresponding field strengths. 

As in the Nielsen- Olesen situation we make a change of variables to clarify the physical 
content of the theory. Thus we set 

<f> = Re i0 (5.6a) 

cr = 5e’ x (5.65) 

= V„ - v„0) (5.6c) 

An = ~{Qn ~ (5.6d) 


In terms of these variables the new lagrangian becomes 

C = V„iZV"iZ + R 2 P^P» - + V m 5V"S + S 2 Q fi Q lt - 

4 g* 4e* 


^-{R 2 - V 2 ) 2 + ( fR 2 - m 2 )S 2 + ^S 4 


A« 

4 


(5.7) 


where G^ v and F^ v are now the field strengths associated with P M and respectively. 
The associated equations of motion are 


V^iZ - iZP^P" + ~-R(R 2 - n 2 ) + fRS 2 = 0 

(5.8a) 

V m V“5 - SQ^Q* + ^5 3 + (/fZ 2 - m 2 )S = 0 

(5.85) 

+ g 2 R 2 P v = 0 

(5.8c) 

+ t 2 S 2 Q v = 0 

(5.8<f) 


In analogy with the Nielsen-Olesen vortex, we solve equations (5.8) for a ‘static’ 
cylindrically symmetric superconducting string, i.e. one with constant current in the z- 
direction, \ oc z. We will write this constant of proportionality as £. This means that we 


19 


ORIGINAL PAGE 13 
OF. POOR QUALITY 


can choose a gauge in which 


R = R(p) , 5 = S(p) , 

P„ = P(p)V ^ = <?(p)V„z 


(5.9) 


in cylindrical polar coordinates{p,r,tf}. The equations of motion thus become 

-R" - — +**?- + ^R(R 2 - rj 2 ) + fRS 2 = 0 

p p 2 2 

-S” ~ — + SQ 2 + % S 3 + ( fR 2 - m 2 )S = 0 

p 2 


(5.10a) 

(5.106) 


-P" + — + g 2 R 2 P = 0 
P 

-Q" - — + e 2 S 2 Q = 0. 

P 


(5.10c) 

(5.10<f) 


We need to impose suitable boundary conditions for a vortex solution. These are 


R —* 0 P -* 1 Q —* C S —> S 0 as r — >0 
R — *7] P — > 0 5 — » 0 as r — ► oo 


(5.11) 


where So is given in terms of the other constants. These have been examined numerically 
in refs [18]. However, we need the analytic expressions for the solution near the origin. 
These can be seen to be 


Pocp , P = l + 0(p 2 ) 

S = S 0 + 0(p 2 ) , Q = C + 0(p 2 ). 


(5.12) 


Equation (5.10b) is satisfied by appropriate choice of the 0(p 2 ) coefficient of Q. Therefore 
we can read off the expression for the gauge fields in this vortex model as 


c, = \(P(p) ~ 1)V„9 

K = |(Q(p) - 


(5.13) 
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As with the Nielsen- Olesen vortex, the gauge fields (5.13) modify the Klein-Gordon 
equation (5.2). The radial equation now becomes 


i 0,p9Mf>) = [(«(/>) - <) ! + (P(p) - 1 ) 2 /p 2 ] *>(/>) 

= 0{p*)v(p) 


(5.14) 


which allows <p(p) —> const, as p — ► 0. Therefore, as with the Nielsen-Olesen vortex, on 
the scale of the core of the string, we do not have suppression of charged particle wave 
functions. 

In order to calculate the baryon current we require the expression of in spherical 
polar coordinates. This is 

A r = ~(<?(r sin B) — 1) cos 6 

e 

Ag = — ~(Q{r sin#) — l)sin0 (5.15) 

e 

A(f, — 0 . 


From (2.9b) we can see that B ** — B £ o , the baryon current for the ordinary (Nielsen- 
Olesen) cosmic string, since the gauge field in (5.15) has no 4> dependence or <f> dependence. 
Therefore the radial kink cannot carry baryon number in this case. But, as with the 
Nielsen-Olesen vortex we will consider an unwinding of topological charge where all three 
pion fields approach the core of the string. 

As before we use the time dependent nucleon Ansatz (4.18). The calculation of the 
baryon current proceeds in a similar fashion to the Nielsen-Olesen case. Since B M = B £ 0 , 
the radial baryon current is given by (4.19) as before. Thus we get the same baryon flux 
(4.21) as with the ordinary cosmic string, 

/ B r d 2 x = -(1 - cos 2 F), (5.16) 

Js* * 

giving the same baryon number flux as before, and hence 


A N b = -1 
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as before. 

Thus superconducting strings catalyse baryon decay. But, since we were forced to 
take a vortex model, i.e. a string with thickness, the process proceeds on the scale of the 
string core. Note once again, the important feature that on intermediate scales (between 
the core radius and the Compton wavelength of the nucleon), the Nielsen- Olesen gauge 
field dominates the behaviour of the ^-field. Thus any suppression of charged particle wave 
functions is the same as for ordinary non-superconducting cosmic strings. 
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6. Discussion and conclusions. 

We have developed the argument of Callan and Witten for monopole catalysis of 
skyrmion decay in such a way that the effects of a topologically non-trivial gauge field are 
highlighted. We then explained the corresponding scenario for cosmic strings. We found 
that a wire model of the string was incompatible with catalysis, but that a vortex model 
admitted a catalysis scenario. This was also shown to be the case for superconducting 
strings. For the superconducting string we first had to derive the vortex solution, since 
previously this had only been found from numerical studies [18]. We did this by solving 
the string equations of motion in the core. With our vortex model we then showed that 
catalysis occured. 

These results support the following heuristic argument for the enhancement factor in 
the case of the monopole cross-section. Notice that the monopole argument was conducted 
exclusively within the approximation of the Dirac monopole; the only place the concept of a 
grand unified monopole occurred was in invoking baryon number non-conserving boundary 
conditions. By contrast, a thick string or vortex model was required in order to get catalysis 
to occur at all in the string picture. Thus in the monopole picture, the only scale we have 
is the skyrmion scale — giving an expected cross-section of the order of a strong interaction 
cross-section. On the other hand, the inescapability of the vortex model in the string case 
suggests that the reaction is occurring on the scale of the string radius, rather than the 
skyrmion radius, thus giving a grand unified cross-section. 

So far our analysis suggests that the monopole catalysis cross-section is proportional 
to m -1 whilst the string and superconducting string cross-sections are proportional to 
M _1 , where m and M are the fermion and grand unified masses respectively. However, 
it is also possible to give the order of magnitude for the constant of proportionality. For 
the monopole we would expect the constant to be 0(1). This is because the process 
proceeds via the radial kink, and there is no suppression of the neutral pion wave function 
in the presence of the monopole. For the case of the string we can find the constant by 
an examination of the Klein-Gordon equation. For distances between m -1 and M' 1 the 
equation for the charged particle wave function is (4.4). Thus, for these distances ~ p. 
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But for p < M~ l the relevant wave equation is (4.15) and ip ~ const, as p — > 0. To match 
solutions at p = M we require that amplitude be proportional to m/M. Thus, we would 
expect this extra suppression factor for the string, giving a cross-section proportional to 
. Similar arguments apply to the superconducting string case. 

These results support the earlier calculations involving a quark/string scattering, that 
is, that there is no enhancement of the baryon decay cross-section for strings. Hence there 
will be no constraints on the cosmic string scenario from catalysis based on later time 
astrophysical processes. However, at earlier times in the evolution of the universe the 
string distribution was more concentrated. It is possible that this could have influenced 
the development of baryon asymmetry. 

We realise these arguments are incomplete. A calculation of the cross-sections in each 
case (string and monopole) is required. This work is in progress. However, as it stands, the 
argument provides an elegant pictorial description of the skyrmion decay process. It shows 
clearly the difference between the monopole and string cases, and also readily obtains the 
superconducting string catalysis picture. 
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Appendix: Calculation of the radial baryon current. 


In this appendix we will calculate the baryon current for the Ansatz nucleon field 

(4.18) 

Us{x,t) = exp[i.P(r,<)i.T] 

= cos F + t sin Fx.r . ( Al) 

We will use spherical polar coordinates, in which the cosmic string gauge field takes the 
form 


A„ = i(P(r sin 0) - 1)V M 0 . (A2) 

In deriving the baryon current 

B „ = {u^d^UU^d^U^d^ (.) 

+ ^ e^d v {A P ^Q{U^d a U + drUU-')) (b), (A3) 

the following expressions will be useful. Recall the identity 

TiTj = 6ij + itijkTk • (A4) 

This implies 

TiTjTk = TiSj k - TjSik + T k Sij + itijkl, • (A5) 

We also need the expression for the derivative of Us, 

Us = iFx.rUs (A6a) 

* * El 

diU N = ix.TF ti U N + ~~~(t* - iix.r) . (A66) 

T 

Consider firstly term (•) in the baryon current (A3). From (A6) we may deduce 

U^doUs = iFx.r (A7a) 

UjfdiUs = ~ [r£jF t i + sin F cos F (6ij - Xi&jj + sin 2 . (A76) 
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For notational convenience, we will write (A7b) as 

U^diU N = iTjMji . 

Then the radial part of (•) is given by 

(.)„ = {y^a v v N v^a p u N v^a,v N ) 


(AS) 


247T 2 

—i 
247T 2 

iF 


4t r 2 r 2 sin# 

F 


ZiMjgMk4>Tr(TiTjTk ) 


2tt 2 


(MggM^ - Mj^) 


F sin 2 F 
-~2^ r 2 

In order to calculate the radial current arising from term (b), we note that 

«-(« A) 

=> TrQri = #< 3 . 


(A 9) 


(A 10 ) 


Thus 


Tr Q(U^Un + U N U j?) = 2 iFiiTrQn 

= 2iFiz 


(ah) 


and 

Tr Q(U^diU N + diU N U = TrQ(2ix.TF ti + 2i sin F cos F(6ij - 

= 2ix%F t i + 2zsinFcosF(#j3 — XiXj)lr (-A12) 

imply 

TV Q(UZ 1 V n + V n U„ 1 ) = 2iF cos 9 ( A13a) 

Tr Q{U^ U N ,g + U N , e U ^ ) = -2 i sin F cos F sin 9 . (A136) 
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Therefore the contribution to the radial baryon current from term (b) becomes 


(b) r = ^2'r UP<T dv[ApTrQ(V- 1 d < rU + d < ,UU- 1 )} 

— ^j-e r °^|5o[— 2i sinFcos FsinflA^] — S^t'Fcos j 


1 F 


47r 2 r 2 sin 


t^|[— cos 2F sin 6(P(r sin 6) — 1)] — 0$[cos 8(P(r sin 6) — 1)]| 


f- cos2F(P - l) + (P - 1) - rP'^4- ) 
l sin 9 ) 


4rr 2 r 2 


{A 14) 


Combining the two expressions (A9) and (A14) gives the total radial baryon current 
Psin 2 P F f „ _ v -v —i cos 2 6 ' 


B r = 


27 r 2 r 2 47r 2 r 2 

F 


{ - cos 2 F(P - 1) + (P - 1) - rP'^-^-) 
l* sin 0 J 


4ir a r 2 


[P(cos 2P - 1) + rP 1 COS j l . 
L sin# J 


(- 415 ) 
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